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The sensitivity to the sterile neutrino mixing at very short baseline 
reactor neutrino experiments is investigated. If the reactor core is relatively 
large as in the case of commercial reactors, then the sensitivity is lost for 
Am 2 > 1 eV 2 due to smearing of the reactor core size. If the reactor core is 
small as in the case of the experimental fast neutron reactor Joyo, the ILL 
research reactor or the Osiris reactor, on the other hand, then sensitivity 
to sin 2 2#i4 can be as good as 0.03 for Am 2 ~ several eV 2 because of its 
small size. 

PACS numbers: 14.60.Pq,25.30.Pt,28.41.-i 

1. Introduction 

Schemes with sterile neutrinos have attracted a lot of attention since 
the LSND group announced the anomaly [H [H [3] which would imply mass 
squared difference of 0(1) eV 2 if it is interpreted as a phenomenon due to 
neutrino oscillation. The standard three flavor scheme has only two inde- 
pendent mass squared differences, i.e., Am^ = Am?) ~ 8 x 10 _5 eV 2 for the 
solar neutrino oscillation, and [Am 2 ,! = Am 2 tm ~ 2.4 x 10 _3 eV 2 for the 
atmospheric neutrino oscillation. To accommodate a neutrino oscillation 
scheme to the LSND anomaly, therefore, the extra state should be intro- 
duced. This extra state should be sterile neutrino, which is singlet with 
respect to the gauge group of the Standard Model, because the number of 
weakly interacting light neutrinos should be three from the LEP data [1] . 

Recently sterile neutrino scenarios are becoming popular again because 
of a few reasons. One is the data of the MiniBooNE experiment which been 
performed to test the LSND anomaly. Although their data on the neutrino 
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mode [5J disfavors the region suggested by LSND, their data on the anti- 
neutrino mode [6] seems to be consistent with that of LSND. The second 
one is the so-called reactor anomaly. The flux of the reactor neutrino was 
recalculated in Ref. [7] recently and it was claimed that the normalization is 
shifted by about +3% on average. This claimed is qualitatively consistent 
with an independent calculation in Ref. [SJ. If their claim on the reactor 
neutrino flux is correct, then neutrino oscillation with Am 2 > leV 2 may be 



concluded from a re-analysis of 19 reactor neutrino results at short base- 
lines [9] • The third one is the so-called gallium anomaly. The data of the 
gallium solar neutrino calibration experiments indicates deficit of v e and it 
may imply neutrino oscillation |10j . 




It has been known that reactor experiments with more than one detector 
has a possibility to measure #13 precisely because some of the systematic 
errors can be canceled by the near- far detector complex [12| [131 E] • 
Three experiments |15|. \T7\ [16] are now either running or expected to start 
soon to measure #13. In the standard three flavor case with (Am^j = 
2.4 x 10 _3 eV 2 , it was shown assuming infinite statistics that the optimized 
baseline lengths Lp and Ljv for the far and near detectors are Lp ~ 1.8km 
and Ln —0km in the rate analysis |18[ 119] . while they are Lp ~10.6km 
and Ltv —8.4km in the spectrum analysis [20]. To justify the assumption 
on negligible statistical errors for L ~10km, unfortunately, one would need 
unrealistically huge detectors, so one is forced to choose the baseline lengths 
which are optimized for the rate analysis for Am 2 = 2.4 x 10 _3 eV 2 . On the 
other hand, if one performs an oscillation experiment to probe Am 2 ~ 0(1) 
eV 2 , it becomes realistic to place the detectors at the baseline lengths which 
are optimized for the spectrum analysis (See Sect. 4 in the published version 
of Ref. [2D]). 

In this talk I would like to discuss the sensitivity of very short line reactor 
experiments to the sterile neutrino mixing for Am 2 ~ 0(1) eV 2 in the so- 
called (3+l)-scheme [ST]. Proposals have been made to test the bound of 
the Bugey reactor experiment [22J on the sterile neutrino mixing angle using 
a reactor [251 1210, an accelerator [281 [27] , and a /3-source [291150] . 



Four-neutrino schemes consist of one extra sterile state and the three 
weakly interacting ones. The schemes are called (3+1)- and (2+2)-schemes, 
depending on whether one or two mass eigenstate(s) are separated from 
the others by the largest mass- squared difference ~ 0(1) eV 2 . The (2+2) 
scheme is excluded by the solar and atmospheric neutrino data[31j, so I 

1 See, e.g., Refs. (the published version), [251 126j for earlier works on search for 
sterile neutrinos at a reactor. 
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will not discuss the (2+2) schemes here. In the (3+1) scheme, the phe- 
nomenology of solar and atmospheric oscillations is approximately the same 
as that of the three flavor framework, so there is no tension between the so- 
lar and atmospheric constraints. However, the (3+1) scheme has a problem 
in accounting for LSND and all other negative results of the short baseline 
experiments. To explain the LSND data while satisfying the constraints 
from other disappearance experiments, the oscillation probabilities of the 
appearance and disappearance channels should satisfy the following rela- 
tion 1521135]: 

1 

4 



sin 2 26> LSND (Am 2 ) < - sin 2 26>B U gey(Am 2 ) • sin 2 26> CDHS w(Am 2 ) (1) 



where 6* LSND (Am 2 ), 6* CDHSW (Am 2 ), 6 , Bugey(Am 2 ) are the value of the effective 
two-flavor mixing angle as a function of the mass squared difference Am 2 in 
the allowed region for LSND {pa — > u e ), the CDHSW experiment |34| (z/„ — > 
v^), and the Bugey experiment [22] {p e — > P e ), respectively. The reason that 
the (3+l)-scheme to explain LSND has been disfavored is because Eq. ([1]) 
is not satisfied for any value of Am 2 , if one adopts the allowed regions in 
Refs. [34J and [22] . If the flux of the reactor neutrino is slightly larger than 
the one used in the Bugey analysis [22J, however, the allowed region becomes 
slightly wider and one has more chance to satisfy Eq. ([T]H. 

I will use the following parametrization for the mixing matrix [ 



U — i?34(#34, 0) i?24(^24, 0) i?23(^23 5 fa) #14 (#14, 0) #13 (#13, 62) #12 (#12, #l) , 

where Rjk(8jk, &i) are the complex rotation matrices in the j/c-plane defined 
as 

[Rj k (9jk, $l)]pq 

= S pq + (cose jk - l)(5jp5 jg + S kp S kq ) + s'm6 jk (e~ lSl 5 jp 5 kq - e l6l 5 jq 5 kp ). 

With this parametrization, for the very short baseline reactor experiments, 
where the average neutrino energy E is approximately 4MeV and the base- 
line length is about 10m, I have \/\m?- k L/AE\ <C 1 (j, k = 1,2,3), so that 
the disappearance probability is given by 

P(v e -+9 e ) = l- sin 2 2^i4 sin 2 f^f^) (2) 

to a good approximation. So the analysis of the (3+l)-scheme is reduced 
to that of a two flavor framework with the oscillation parameters (Am^, 
sin 2 26»i4 ). 



2 Although the situation of the (3+l)-scheme is improved slightly after Refs. [5], 
the improvement is not sufficient enough to have a satisfactory fit to all the data, 
according to Ref. [35] . 
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3. Sensitivity to sin 2 2#i4 by a spectral analysis 

Throughout my talk I discuss the case with a single reactor and two de- 
tectors. I assume here that the near and far detectors are identical and they 
have the same sizes of systematic errors. The conditions of the detectors 
are assumed to be the same as those of the Bugey experiment, i.e., liquid 
scintillation detector of volume 600 liters with the detection efficiency which 
yields about 90,000 events at L=15m from a reactor of a power 2.8GW after 
running for 1800 hours. 

To evaluate the sensitivity to sin 2 29\4, let us introduce the following \ 2 
which was adopted in Ref. [20] (See Ref. [21] for details): 



ri 
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A=N,F 



(3) 



X 2 stands for a quantity which expresses how much deviation we have be- 
tween the numbers of events with and without oscillations, compared with 
the experimental errors. In Eq.([3j), mf is the number of events to be mea- 
sured at the near (A = N) and far (A = F) for the i-th. energy bin with 
the neutrino oscillation, and tf is the theoretical prediction without the 
oscillation, (erf-) 2 is the uncorrelated error which consists of the statistical 

plus uncorrelated bin-to-bin systematic error: (tferf) 2 = tf + (tf&fb) , 

where a A h is the uncorrelated bin-to-bin systematic error. a A (A = N, F) 
is a variable which introduces the detector-specific uncertainties <7<3B of the 
near and far detectors, (i = 1, • • • ,n) is a variable for an uncertainty 
o"Db of the theoretical prediction for each energy bin which is uncorrelated 
between different energy bins|f| af^ (A = N, F) is a variable which intro- 
duces an energy calibration uncertainty <r ca i and comes in the theoretical 
prediction in the form of (1 + af^E instead of the observed energy E. vf is 
the deviation divided by the expected number of events from the theoretical 
prediction tf due to the energy calibration uncertainty. Here I take the fol- 
lowing reference values for the systematic errors: a^b = 0.5%, <7dB = 0.5%, 
a Bh = 2%, ctdb = 3%, a cal = 0.6%. 



3 The first suffix of a stands for the property for the systematic error with respect 
to the detectors while the second is with respect to bins, and capital (small) letter 
stands for a correlated (uncorrelated) systematic error. 
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3.1. Commercial reactors 

First of all, I will consider a commercial reactor whose thermal power is 
2.8GW and I will assume that the dimension of its core is 4m in diameter 
and 4m in height. 

X 2 in Eq. ([3]) is computed numerically in the case of Am^ = leV 2 as 
a function of the baseline lengths Ln and Lp of the two detectors, and 
the baseline lengths Ln and Lp are varied to optimize the sensitivity to 
sin 2 2#i4. It is found that the set (L]y,Lp) ~ (17m, 23m) gives the optimum. 
In contrast to the rate analysis, in which the optimized baseline length of 
the near detector is LAr=0m to avoid oscillations, the spectrum analysis 
with (L]y,Lp) = (17m, 23m) looks at the difference between the maximum 
and minimum of the spectrum shape with neutrino oscillations at Ln and 
Lp mainly for the energy region E u ~ 4MeV where the number of events are 
expected to be the largest (See the upper panel in Fig.[TJ). Unlike the case 
of infinite statistics |20| , the statistical errors are important in the present 
setup of the detectors, and longer baseline lengths are disfavored. 

The sensitivity to sin 2 29^ in the case of the baseline lengths {L^,Lp) = 
(17m, 23m) is shown in Fig. [2] as a function of Am]] (the line referred to 
as "Commercial"). The region suggested by combination of the reactor 
and gallium anomalies and the MiniBooNE data is also given in Fig. [2] for 
comparison. For Am^ >2eV 2 , the sensitivity is no better than 0.1, which 
is basically the result of the rate analysis. The sensitivity in the case of a 
hypothetical point-like reactor, where all the conditions for the detectors 
are the same, is also given in Fig. [2] for comparison (the line referred to 
as "Point-like"). Fig. [5] indicates that the sensitivity would be as good as 
several x 10 -2 for a few eV 2 , if the core were point-like. So we can conclude 
that we have poor sensitivity for Am 2 ] > 2eV 2 because of the smearing 
effect of the finite core size of the reactor. 

3.2. Research reactors 

In the previous subsection, we have seen that the sensitivity to sin 2 2#i4 
is lost because of the smearing effect of finite core size. Next, I would like 
to discuss three research reactors, Jovo [37] with MK-III upgrade [38], the 
ILL research reactor [39] , and the Osiris research reactor [30] . They all have 
a relatively small size and a relatively large thermal power. 

Joyo is an experimental fast breeder reactor and the dimension of its core 
is 0.8m in diameter and 0.5m in height, and its thermal power is 140MW. 
The ILL (Osiris) research reactor is a thermal neutron reactor with high 
enrichment uranium 235 U, and the dimension of its core is 0.4m in diameter 
and 0.8m in height (0.57mx0.57mx0.6m) and its thermal power is 58MW 
(70MW), respectively. 
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Fig. 1. The energy spectrum with neutrino oscillations at the two different detec- 
tors and the one without oscillations. The optimized baseline lengths give maxi- 
mum difference in the distortions in the energy spectrum. The upper panel: the 
case of a commercial reactor. The lower panel: the case of a research reactor. 

Again \ 2 m Eq. © is computed numerically in each case, and it is 
optimized with respect to Ln and Lp. The optimum set of the baseline 
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Joyo with L F =4m, L N =3m 
Joyo with L F =8m, L N =4m 
ILL with L F =8m, L N =4m 
Osiris with L F =8m, L N =4m 
Commercial with L F =23m, L N =17m 
Point-like with L F =23m, L N =17m 

Reactor anomaly 90%CL 




sin 2 26 14 

Fig. 2. The sensitivity to sin 2 29\^ of each reactor with the two detectors at its 
optimum baseline lengths. Also shown as a shaded area is the region given in 
Ref. [9] from the combination of the reactor neutrino experiments, Gallex and Sage 
calibration sources experiments, the MiniBooNE reanalysis of Ref. [10] . and the 
ILL-energy spectrum distortion. 
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lengths turns out to be {L^,Lp) ~ (4m, 8m) for Am^ = leV 2 for all the 
three cases. The lower panel in Fig. Q] shows the spectrum distortion in the 
case of L=4m, 8m. 

The sensitivity to sin 2 2#i4 is shown in Fig. [2] as a function of Am^ in 
the case of the sets of the baseline lengths (L^,Lp) = (4m, 8m) for the 
three cases and {L^,Lp) = (3m, 4m) for Joyo. From Fig.[2]it is clear that 
the sensitivity of an experiment with a small core reactor is better that that 
with a commercial reactor for 2eV 2 <Am1i <10eV 2 . 

4. Discussion and Conclusion 

In the framework of the (3+l)-scheme, the sensitivity to sin 2 2#i4 of very 
short baseline reactor oscillation experiments was studied by a spectrum 
analysis. The assumptions are that one has two detectors whose size and 
efficiency are exactly the same as those used at the Bugey experiment and 
X 2 is optimized with respect to the positions of the two detectors. 

In the case of a commercial reactor, by putting the detectors at Ln = 
17m and Lp = 23m, one obtains the sensitivity as good as several xl0~ 2 
for Am 2 ! <leV 2 , but the sensitivity is lost above leV 2 due to the smearing 
of the finite core size. 

In the case of a research reactor with a small core (such as Joyo, ILL, 
Osiris), on the other hand, one obtains the sensitivity as good as a several 
xlO -2 for leV 2 <Am?, <10eV 2 if the detectors are located at Ln = 4m 
and Lp = 8m. 

In all the cases discussed above with the Bugey-like detector setup, the 
statistical errors are dominant. The reason that the case of the research 
reactors (Joyo, ILL, Osiris) is competitive despite its small power is because 
the total numbers of events at L ~ several meters are comparable to those 
of the case with a commercial reactor at L ~ a few x 10 meters. 

To turn this idea into reality, there are two experimental challenges. One 
is to put detectors at a location very near to a research reactor. The other 
one is to avoid potentially huge backgrounds from the reactorjf] 

Nevertheless, since the best fit point (Am^, sin 2 29 u ) ~ (2eV 2 , 0.1) 
obtained in Ref. [S] lies within the excluded region in Fig. [21 the experiment 
at these research reactors offers a promising possibility. 
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